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Abstract. In this work we performed density functional theory calculations for 
a twisted bilayer graphene (BLG). Several conmensurable rotation angles were 
analyzed and for each one a constant height mode STM image was obtained. These 
STM images, calculated under the Tersoff-Hamman theory, reproduce the main 
features experimentally observed, paticularly superstructures and giant corrugations. 
In this way we confirm that STM characterization of twisted BLG can produce 
superstructures whose tunneling current intensity maxima occur over regions with 
AA stacking. Additionally we give new evidence in favour of an electronic origin for 
the superstructures instead another physical grounds. 
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1. Introduction 

The increasing interest for nanoscience and nanotechnology has evidenced the 
importance of the scanning tunneling microscope (STM) to explore the nature at atomic 
scale. In this context graphene, a novel layered material with potential technological 
applications [1] , has been intensively characterized by STM techniques and has revealed 
its surprising behavior: edge states at the border [2j [3], superstructures [U El |6] 
and Van Hove Singularities [7]. The graphene particularities were implied before 
by STM characterization of the (0001) surface of highly oriented pyrolitic graphite 
(HOPG): difficulties to obtain STM complete atomic resolution as well the observation 
of superstructures are examples of this. 

Extending the preceding statements, let us precise that graphene is a stable 
bidimensional structure formed by carbon atoms ordered in a hexagonal lattice (the first 
neighbor distance is 1.42 A). In crystalline graphite (HOPG) such layers are stacked on 
each other in a AB sequence, and therefore it is possible to identify two type of atomic 
sites: the a-type sites, corresponding to atoms with neighbors directly above and below 
in adjacent layers; and the /9-type sites, which do not have direct neighbors in the 
adjacent layers. An important consequence of the crystalline structure of HOPG is that 
near the Fermi level the bulk presents a Local Density of States (LDOS) which is larger 
for atoms on /3-type sites [8]. The theory of STM image formation [9j[T0] considers that 
the accessible states by STM are those in the energy window [E Fermi — eVu aS) E Fermi ]. 
Consequently in HOPG only those atoms on /3-type sites appears visible for a STM 
operating at low bias voltages and they form the triangular lattice usually reported on 
experiments. However the LDOS of a and (3 sites tends to be equal by increasing the 
bias voltages and the persistence of the triangular lattice at higher bias voltages has 
been explained by resolution losses introduced for the STM tip size [TT] . 

Following the preceding ideas, the observation of superstructures during STM 
experiments can be explained by considering a relative rotation among the constituent 
graphene layers. This is because rotations generate Moire Patterns where it is possible 
to identify regions with a high concentration of /3-type sites, called g — (3 sites, and 
regions with a high concentration of a-type sites, called g — a or AA stacked regions 
[see Figure (Ha)]. In both cases g stands for giant [12]. In this frame the g — (3 regions 
would appear brilliant in STM images due to the large LDOS of atoms on (3 sites 
[T2] . However, Rong and Kuiper [13], on base of first principles DOS calculations [H], 
proposed that brilliant zones occur over the g—ct regions. Aditionally a third explanation 
for the superstructures considers that more important than the electronic effects induced 
by rotation between weakly interactive graphene layers, could be surface deformations 
introduced by STM tip while it scans the surface (see Ref. [15] and references therein). 
This has been a long standing controversy which has attracted renewed interest over the 
STM image formation of graphite surface [HJ [HI El EEB] and over the physical origin of 
superstructures [TBI [i~9l [20l |2T] . 

In this context, an important tool to analyze the effect of rotation among 
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graphene layers is the calculation of STM images and the corresponding contrast with 
experimental data. Following this objective we have focused in the task of generate STM 
images for different rotation angles between two graphene layers. We have organized this 
paper as follows: details of the images calculation are given in next section; theoretical 
results are summarized in section [3] and finally the main conclusions are presented. 

2. Calculation Method 

The system under study corresponds to BLG presenting a relative rotation angle between 
its layers (twisted bilayer graphene). Such rotation occurs around the stacking direction 
and despite for any angle one can identify a superstructure (or moire pattern) only for 
very particular angles the misoriented layers are in commensuration [2TI 122] [23] . A 
commensurable unit supercell requieres a rotation from a vector V\ = max + na,2 to 
V 2 = na± + ma 2 , where a\ = (\/3, — l)ao/2 and a 2 = (V%, l)oo/2 are the graphene 
basis vectors; m and n are integers and a = 2.46 A is the lattice constant. The 
commensurable rotation angle is defined by 

m 2 + 4mn + n 2 

COSfc>=— — 5 57-, (1) 

2(777/ + 77777 + 77^) 

and the unit super cell vectors correspond to: t\ — V 2 — na\ + 777a 2 and t 2 = 
—di + (m + n)a 2 [21] • The unit supercell contains iV = 4(m 2 + mn + n 2 ) atoms and 
its periodicity results D = a o /[2sin(0/2)] [121 [15]. Thus, commensurable unit supercells 
are complety defined by the index m and 77 and those selected to this study appear in 
Table [1] with their corresponding number of atoms N, rotation angle between layers 9 
and periodicity D. For example, the unit supercell for twisted bilayer graphene with 
the angle 7.3° is shown in Figure 1(b). 

Table 1. Selected commensurable unit supercells and their corresponding number of 
constituent atoms (N), rotation angle (8) and periodicity (D). The last column shows 
periodicities experimentally reported. 



Ill 


11 


N 


8 


D (nm) 


D Exp. 


(nm) 


2 


1 


28 


21.8° 


0.65 


0.65 Ref. [H] 




3 


2 


70 


13.2° 


1.05 


1.06 Ref. [15] 


0.95 Ref. |25] 


4 


3 


148 


9.4° 


1.50 


1.50 Ref. [H] 


1.50 Ref. [1H] 


• r j 


4 


244 


7.3° 


1.91 


1.76 Ref. [IS] 


1.71 Ref. [3D] 


(i 


5 


364 


6.0° 


2.30 


2.20 Ref. [31] 


2.40 Ref. [32] 



Our DFT calculations were carried out using the SIESTA ab initio package 
[33] which employ norm- conserving pseudopotentials and localized atomic orbitals as 
basis set. Double-^ plus polarization functions were used under the local density 
approximation [M]. All structures were fully relaxed until the atomic forces are smaller 
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than 0.02 eV/A. We consider supercells with periodic boundary conditions for the plane 
layers, while in the perpendicular direction to the layers plane the separation between 
slabs is fixed to 14 A. The Brillouin zone sampling was performed using a Monkhorst- 
Pack mesh of 10 x 10 x 1. 

The images were obtained using the code STM 1.0.1 (included in the SIESTA 
package). This code uses the wavefunctions generated by SIESTA and computed on a 
reference plane and extrapolates the value of these wavefunctions into vacuum. Such 
reference plane must be sufficiently close to the surface so that the charge density is 
large and well described. The STM data is generated under the Tersoff-Hamann theory 
[9], while data visulization was possible using the WSxM 5.0 freeware [35] . Gaussian 
smooth was performed to obtain the final STM image. 

3. Results 

Electronic properties and constant height mode STM images were calculated for the five 
commensurable supercells presented in Table [TJ Figure [2] shows band structure and total 
density of states (DOS) for three of these twisted bilayer graphene. Clearly, a decrease 
of the rotation angle 9 diminishes the energy difference between the VHS near the Fermi 
Level. This effect was reported in previous works [TJ [21j [36] and has its origin in the 
renormalization of Dirac electrons velocities, which is observed in the band structures 
as a decreasing in the slope of the Dirac cone. Additionally the VHSs show the well 
known asymmetry between electrons and holes, which is observed as a difference in the 
height of the VHSs. These theoretical predictions have been recently corroborated by 
STM/STS experimental measurements on twisted graphene layers with 9 = 1.79° [7] 
and ensure us that our calculations reproduce the main electronic characteristics of the 
twisted bilayer graphene. 

As the occupied states contributing to the STM image are in the energy window 
[Epermi — e\ Vbi as \ , Epermi] [3 EH] , the STM image bias voltage dependence was studied 
by considering several values in the range 0.05 V < Vu a s < 3.0 V. As example, 
Figure [3] shows calculated STM images and their corresponding line profiles for the 
case Vbi as — 1-0 V and the tip-to-surface distance oIts — 1-0 A. These calculated STM 
images show bright zones indicating tunnel current maxima, which form triangular 
structures with periodicities much larger than the graphite lattice constant. These 
superlattice constants correspond to the periodicities D shown in Table [H Another 
important result derived from these calculations is that the brilliant zones are over the 
AA stacked regions in accordance with the proposal of other authors: Rong and Kuiper 
[13] and Campanera et al. [19], who based their conclusions in DOS calculations for 
different stacking sequence of graphite |14j ; and Trambly de Laissardiere et al. |21j . 
who combined tight-binding with ab-initio calculations to study the LDOS on twisted 
bilayer graphene. 

Line profiles were performed along the segment indicated in each STM image of 
Figure [31 Although the calculated images correspond to constant height mode, the line 
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profiles might be compared with those experimentally obtained in constant current mode 
assuming that depressions in the charge density lead to proportional approximations of 
the tip-to-the-surface in order to keep the tunnel current constant. So, our results 
show two important features revealed also during the STM experiments: a roughnes 
presenting the graphite periodicity and a large oscilation amplitude presenting the 
superlattice periodicity [T2 | IT3 | IT5] . For the first case, the roughnes reflects the electronic 
charge accumulation at the atomic sites, which is detected by STM tip. For the second 
case, the oscillation amplitude is much more larger than the roughness due to the atomic 
sites and it has been named "giant" or "superlattice corrugation" [15]. Clearly our 
calculations suggests that the physical origin of the superlattice corrugation is mainly an 
electronic effect and consequently the wide gamma of alternative explanations (see [15] 
and references therein) seems to be less important at least for STM images of BLG. Also 
we found, coinciding with the experimental observations on graphite surfaces [131 US] , 
that bias voltage has a strong influence in the Superlattice Corrugation Amplitude 
(SCA). Figure H] shows SCA obtained from calculated STM for the rotation angles 
under study and with different bias voltages. The SCA corresponding to 9 = 21.8° 
is the smaller one and presents a current maximum at ~ 1.0 V. Such SCA maximum 
seems to be a consequence of the DOS which, in the case 9 = 21.8°, shows the first 
VHS at Rj —1.0 V. Of course the main contribution to this VHS comes from the Local 
Density of States (LDOS) over the AA stacked regions, but the decrease of SCA for 
Vbias > 1-0 V is evidence that from that value other regions start to contribute to the 
DOS. For lower rotation angles, as i.e. 9.4° and 7.3°, the first VHS occurs closer to the 
Fermi Energy (see Fig. [2]) and its main contribution comes from the AA stacked regions. 
However when the bias voltage increases over 2.0 V the LDOS from other regions start 
to contribute equally and SCA starts to diminish. Such behavior for the LDOS was also 
experimentally revealed by STS on the bright and dark regions of STM data [3 124"] . 

On the other hand the STM dependence on the tip-to-sample distance (drs) wa s 
also studied by considering different values in the range 1.0 A < g^s < 2.0 A. Thus, the 
present STM image generation methodology shows that the tip-to-surface distance has 
no consequence in the calculated STM images and current intensity maxima remain over 
the regions with large concentration of a-sites (AA regions). This result indicates that 
the STM image generation methodology is well defined in the dxs range considered, but 
does no give definitive conclusions about tip-to-sample distance dependences as those 
theoretically predicted [20] and experimentally observed [6] for trilayer graphene case. 

4. Conclusions 

From ab-initio methods we have calculated constant heigth mode STM images for 
twisted graphene bilayer presenting five commensurable rotation angles. The calculated 
images show superstructures as those revealed during STM characterization of graphite 
surface and reproduce their main features. In particular the line profiles, obtained from 
the calculated images, show two important features experimentally observed: a large 
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oscillation amplitude (giant corrugation) with the periodicity of the superstructures; 
and a roughness with the graphene lattice periodicity. Thus, these results are showing 
that the superstructures have a strong electronic origin and, therefore, mechanical tip- 
sample interaction, or another physical grounds, are not always relevant to explain the 
phenomena. 

Three additional important conclusions derived from this study must be remarked: 
first, despite the weak interlayer interaction, a relative rotation between them can 
modify considerably the surface electronic structure; second, results show clearly that, 
for twisted bilayer graphene, the current intensity maxima appear over AA stacked 
regions; and third, the bias voltage has a strong influence in the superlattice corrugation 
amplitude. 
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Figure 1. (a) Scheme view of twisted bilayer graphene for 7.3° as rotation angle, (b) 
Magnification of the unit supercell which expands the superstructure. Green sticks are 




Figure 2. Bands structure and total density of states for different twisted bilayer 
of graphene. Left panels show the bands structures and right panels show the total 
density of states. 
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Figure 3. Constant height mode STM images calculated for the five commensurable 
unit supercells presented in Table [T] The line profiles (along the line indicated on the 
corresponding STM image) appears at right of each one. 
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Figure 4. Supcrlattice Corrugation Amplitude on calculated STM images as a 
function of bias voltages. The rotation angles for the twisted BLG are indicated in the 
figure. 



